Wannier functions and exchange integrals: The example of LiCu 2 02 
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Starting from a single band Hubbard model in the Wannier function basis, we revisit the problem 
of the ligand contribution to exchange and derive explicit formulae for the exchange integrals in 
metal oxide compounds in terms of atomic parameters that can be calculated with constrained LDA 
and LDA+U. The analysis is applied to the investigation of the isotropic exchange interactions of 
LiCu2 02, a compound where the Cu-O-Cu angle of the dominant exchange path is close to 90°. 
Our results show that the magnetic moments are localized in Wannier orbitals which have strong 
contribution from oxygen atomic orbitals, leading to exchange integrals that considerably differ 
from the estimates based on kinetic exchange only. Using LSDA+U approach, we also perform a 
direct ab-initio determination of the exchange integrals LiCu2 02- The results agree well with those 
obtained from the Wannier function approach, a clear indication that this modelization captures 
the essential physics of exchange. A comparison with experimental results is also included, with the 
conclusion that a very precise determination of the Wannier function is crucial to reach quantitative 
estimates. 

PACS numbers: 73.22.-f, 75.10.Hk 



I. INTRODUCTION 

In order to explain the homeopolar chemical bond in 
the hydrogen molecule, Heitler and London! have derived 
the direct and exchange Coulomb integrals which can be 
written in the following form: 



and 



u 



Coulomb 



dxdx' 



(1) 



\x — x'\ 



dxdx', (2) 



where i and j are site indexes, and <pi(x) is a wave 
function centered at the lattice site i. Taking into ac- 
count the exchange Coulomb integral allows one to ob- 
tain bonding (E = U — jCouiomb^ ag we n as antibonding 

(E' = U+ jc° ut ° mb ) states of H 2 . Heitler and London 
have shown that the properties of the hydrogen molecule 
can be described correctly using these combinations of 
wave functions. 

Then, in 1928, Heisenberg 2 has used the exchange 
Coulomb integral jC° ulomb to explain ferromagnetism. 
Heisenberg has supposed that this exchange Coulomb 
integral corresponds to the exchange coupling in a spin 
model defined by the Hamiltonian: 



H — Jij Si ■ Sj 



(3) 



and that it is the main source of ferromagnetism in 3d 
metal compounds. However, if <f>i{x) in Eq.(l) is the 
atomic wave function centered on the ith atom, then the 
overlap between wave functions on neighboring atoms is 



negligibly small. Therefore the exchange interaction de- 
rived through Eq.(2) can be neglected. 

In 1959, P.W. Anderson 3 - has suggested a new type of 
exchange interaction process based on hopping (kinetic 
exchange interaction). In the context of the Hubbard 
modelii 



ij a tcr a jv 



(4) 



this exchange interaction can be expressed as J^ m = -jf- . 
This estimation of exchange interaction through hopping 
integrals is convenient and now widely used in the liter- 
ature. However, if is small enough, other sources of 
exchange coupling become important. Indeed, the proper 
way to discuss exchange is to consider the basis of Wan- 
nier functions Wi (where i is the composite index of band 
and site) proposed by Wannier— in 1937 and defined as 
the Fourier transforms of a certain linear combination of 
Bloch functions ip n k (n is the band index and k is the 
wave- vector in reciprocal space). In contrast to atomic 
wave functions (f>i(x), which are localized on one atom, 
the Wannier functions Wi(x) are more extended in space 
and can be expressed through linear combination atomic 
wave functions, Wi(x) — Y^jT a j^j( x ~ -0 (where a) 
is the contribution of the jth atomic orbital to Wannier 
function Wi (x) and T is a translation vector) . The Wan- 
nier functions are the most localized ones within the sub- 
space of low-energy excitations, which facilitates a direct 
physical interpretation consistent with interacting local- 
ized spins. Therefore, to use Wannier states instead of 
atomic sets is physically motivated. In particular, as we 
shall see, jy° ulomb defined in the Wannier basis plays a 
crucial role in the description of exchange interactions be- 
tween magnetic moments in the case of nearly 90° metal- 
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oxygen-metal bonds. 

These questions have already been discussed in sev- 
eral context in the literature. For instance, the authors 
of Ref. have discussed two alternative ways, natural 
or orthogonalized magnetic orbitals, to describe the ex- 
change interactions. They have concluded that the or- 
thogonalized magnetic orbital approach clearly leads to 
simpler calculations and thus may be more appropriate 
for quantitative computations. 

The competition between kinetic (</^ m ) and potential 
(jCj oulornb ) contributions to the total exchange interac- 
tion has been considered before in many works. In Ref. 
7, the authors have performed model calculations of a 
CU2O6L14 cluster for different Cu-O-Cu angles. Their 
computational experiments have shown that the nearest- 
neighbor interaction reaches a maximum around 97° and 
remains ferromagnetic up to angles as large as 104°. 
They have also concluded that the simple superexchangc 
relation cannot be applied to Li2Cu02- In Ref. H, it 
was shown that unusual insulating ferromagnetism in 
La4Ba2Cu20io can be explained by intersite ferromag- 
netic "direct exchange" (in our notation jf: j oulomb ). The 
authors have concluded that the latter process occurs 
mainly at La and O sites and overwhelms the AF su- 
perexchange j£? n . The value of jg ouiomb was calculated 
through a direct integration over the wave functions. 

In the present paper, we revisit this issue in the context 
of constrained LDA and LSDA+U ab-initio approaches. 
The problem we want to address can be formulated as 
follows. In favourable situations, the exchange inte- 
grals calculated using LSDA+U agree quite well with 

the standard superexchange expression jf™ 1 — if the 
parameters Uj and U are themselves determined from 
constrained LDA and LSDA+U. Such an interpretation 
of the LSDA+U results is a well accepted criterion to 
test the reliability of the result. But the standard su- 
perexchange expression j£? n = often fails quantita- 
tively, and even sometimes qualitatively, to reproduce the 
LSDA+U results. The main goal of this paper is to pro- 
vide for such cases a generalization of the superexchangc 
expression which is entirely expressed in terms of param- 
eters that can be determined by constrained LDA and 
LSAD+U, and which can be used as an interpretation 
and an independent check of the LSDA+U results. 

To this end, we start from the standard Hamiltonian 
in Wannier function basis for nearly 90° metal-oxygen- 
metal bond. We show how the different Coulomb in- 
teraction terms of this Hamiltonian are related to pa- 
rameters defined in the atomic basis set, with empha- 
sis on the intraatomic exchange interaction of oxygen, 
Jp , to which jCouiomb j g p r0 p 0r tional when neighbor- 
ing Wannier orbitals overlap on the oxygen atoms. We 
then present a simple expression for the exchange inter- 
actions between magnetic moments in the system. The 
parameters that enter this expression can themselves be 
estimated through LDA and LSDA+U calculations. This 
formalism is applied to the investigation of the magnetic 



properties of LiCu202, and the results are compared to 
those of first-principle LSDA+U approach. 

The paper is organized as follows. In Section II, we dis- 
cuss the Hubbard model in Wannier function basis. In 
Section III A, we shortly describe the crystal structure 
of LiCu2C>2 and present the results of LDA calculation. 
In Section III B, we present the results of LSDA+U cal- 
culations and discuss the exchange interactions obtained 
between different pairs of magnetic moments in LiCu202- 
In Section IV, we briefly summarize our results. 

II. HUBBARD MODEL IN WANNIER 
FUNCTION BASIS 

The general Hamiltonian in Wannier function basis 
Wi(x) can be written in the following form: 4 

H = Y t ij a tr a i<r + 2 Y (ij\U\kl)afaf a ,ai a >a k(T ,(5) 

i,j.(T ijkl.aa' 

where 

Uj = J W*{x)V 2 W {x)dx 

and 

r-nrnn f {x)W k {x)W* {x')Wi{x') , 

yij\U\kl) = I ; — dxdx . 

J \x-x'\ 

We analyze the complex Hamiltonian of Eq. (5) in the 
context of a 3-site model (see Fig.l) with nearly 90° 
bonds and define two Wannier orbitals W\ — a<pi + (3<p Pl 
and Wi — acj>2 + l3(p P2 , which are constructed from 
the atomic wave functions: (pi, <p2, <p pi and <p P2 (see 
Fig.l). The oxygen orbitals <p Pl and <p P2 can be expressed 




FIG. 1: (Color online) Schematic representation of a nearly 
90° bond between 3d atoms through oxygen. 

through the angle of metal-oxygen-metal bond 9 in the 
following form: 

Q — go - 90° 
4>pi = cos( )(P Py + sin( )(P Px (6) 

and 

9 - 90° 9 - 90° 
<Pp 2 = cos ( 2 ^ +Sin ( 2 ^ Py ' ^ 
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Let us first analyze the hopping term in Eq.(5). It is 
easy to show that 



tia= / W?(x)VHV. 2 (x)dx = [3 z J (P* pi (x)V 2 (t> P2 (x)dx.(8) 

If 4> Pl and 4> P2 are eigenfunctions of the Hamiltonian of 
Eq.(5), then using the orthogonality condition, one can 
obtain the following expression for the hopping term: 



ti2 ~ /3 2 cos#. 



(9) 



Clearly, if 6 = 90°, the hopping integral vanishes. 

In the analysis of the Coulomb term of Eq.(5) that fol- 
lows, we consider only density-density terms and on-site 
exchange integrals. On-site Coulomb interaction, inter- 
site Coulomb interaction and intersite exchange interac- 
tion in Wannier function basis are expressed through the 
corresponding parameters in the atomic basis set: 



(ii\U\ii) = a 4 U d + 2a 2 (3 2 V pd + (3% 



(ij\U\ij) = a 4 V dd + 2a 2 (3 2 V pd + p 4 U p , 



(10) 



(11) 



(12) 



Ud = / m^m^i dxdxl (13) 



where 



\x — x'\ 

is the on-site Coulomb interaction of 3d atom, 

= f ^(x)^(x)^(xO^(xQ^ 
J |ar — ar'| 

is the on-site Coulomb interaction of oxygen atom, 

V vd = J ^—^\ dxdx ( 15 ) 



which is the so-called correlated hybridization, could sig- 
nificantly change the parameters in the effective single 
band model for transition metal oxides. However we do 
not consider it here and leave that point for further in- 
vestigation. 

Finally one can write the following Hamiltonian in 
Wannier function basis 



{a 4 U d + 2a 2 (3 2 V pd + f3 4 U p ) ^ 

ia 

(a 4 V dd + 2a 2 f3 2 V pd + (3 4 U p ) ^ 

ijcrcj' 



4 tH 



/3 4 J, 



Y,(n ia n j(T +2S?S? + 2SVS%), (19) 



where = a^a J( j is the particle number operator, while 
Sf = |(a^fli| +a+a, T ) and Sf = ^( a t^ a il ~ a ti a il) are 
components of the spin operator. One can reduce Eq.(19) 
to the following form: 



Veff V- 



2 / j %a 3® 



^jfJ^Si-Sj, ( 2 °) 



where the effective on-site Coulomb interaction is given 
by U e ff = a 4 U d +2a 2 f3 2 V pd +/3 4 U p and the effective inter- 
site Coulomb interaction by V e ff = a 4 V dd + 2a 2 2 V pd + 

f3 4 U p y~. It is easy to show that the Heisenberg 

model which corresponds to this electronic Hamiltonian 
has the following form: 



C oulomb 



)Si ■ Sj 



2t 2 



eff ~ V eff 



is the Coulomb interaction between 3d atom and oxygen, where U e tt — V e ff = a> 4 {U d — V dd ) + 



V, 



\x — x'\ 

is the Coulomb interaction between 3d atoms, 



dxdx' (16) 



\x — x'\ 



dxdx' (17) 



is the intraatomic exchange interaction of oxygen atom. 
It was shown in Ref . |9| that the following term: 



(ii\U\ij) = a^ I M X M'W»W*»W 



\x — x'\ 



dxdx' 



-P 4 



<f>* Pl {x)<t>Pi{xW Pl (x')(j) P2 {x') , 



\x — x'\ 



dxdx' (18) 



In the case 

of a nearly 90° metal-oxygen-metal bond (Fig.l), there 
is an additional ferromagnetic contribution to the total 
exchange interaction. The origin of this term is Hund's 
rule exchange interaction on the oxygen atom. As we 
show below, the value of Jy. oulomb is not negligible: it 
can fully compensate the kinetic contribution, so that 
the total exchange interaction becomes ferromagnetic. 

If (3=0, then the Hamiltonian of Eq.(20) is the simple 
Hubbard model. The Coulomb parameters of the gen- 
eral Hamiltonian of Eq.(20) can be calculated through 
constrained LDA calculations or direct integration over 
the wave functions. In this paper we use the constrained 
LDA calculation approach^ which has given reasonable 
results for a number of compounds. We apply the analy- 
sis of this section to the investigation of exchange inter- 
actions in LiCu2 02- 
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III. RESULTS 

A. LiCu 2 2 : LDA CALCULATION 

LiCu202 (Fig. 2) is a good example to demonstrate the 
competition between Coulomb jC° ulomb anc ] kinetic J^ m 
exchange interactions. Indeed, the angle of Cu-O-Cu 
bond in this compound is equal to 94°. This value of 
the angle corresponds to the middle of the ferromagnetic 
range,' and the nearest-neighbor hopping is considerably 
suppressed (see Eq.(9)). 

FIG. 2: (Color online) Crystal structure of LiCu202. Green, 
red, blue and yellow spheres are Cu 2+ , Cu + , O and Li ions, 
respectively. 

D. A. Zatsepin et alM have performed the first ab ini- 
tio calculations of the electronic structure of LiCu2 02 
in terms of LSDA and LSDA+U approximations. They 
have concluded that, in the ground state, LiCii202 is an 
insulator with ferromagnetic ordering. 

A. A. Gippius et al. 12 have performed magnetic reso- 
nance measurements of LiCu2 02 in the paramagnetic and 
magnetically ordered states. They have also performed 
full potential LDA calculations. The band structure near 
Fermi level has been fitted by an extended tight-binding 
model. Using the resulting hopping parameters, the au- 
thors have estimated the values of the exchange interac- 
tions Jij within a single band Hubbard model. 

Experimental investigations of the magnetic ordering 
of LiCu202 by neutron scattering have been performed 
by T. Masuda et aZi 13 i 14 The authors have proposed dif- 
ferent sets of exchange constants obtained by fitting the 
calculated spin wave dispersion relation to the experi- 
mental curves. None of the sets of parameters is in good 
agreement with the results of Refs. [l2|. One can con- 
clude that, at present, there is no consistent theoretical 
and experimental description of the magnetic properties 
of LiCu2 02- 

In the present paper, the electronic structure calcula- 
tion of LiCu2 02 was performed using the Tight Binding 



Linear-Muffm-Tin-Orbital Atomic Sphere Approxima- 
tion (TB-LMTO-ASA) method in terms of the conven- 
tional local-density approximation^ and crystal struc- 
ture data from Ref. [l6|- The band structure of LiCu202 
obtained from LDA calculations is presented in Figj3] 
There are four bands near the Fermi level which are well 
separated from others. These bands are in good agree- 




FIG. 

eV). 



3: Band structure of LiCu 2 02 near the Fermi level (0 



ment with those presented in Ref. 

The partial density of states of LiCu202 obtained from 
LDA calculations is presented in Figj4] Copper 3d states 
of x 2 — y 2 symmetry are strongly hybridized with oxygen 
2p states. Therefore it is more natural to use the Wan- 
nier function basis rather than atomic orbitals to describe 
the hybridization processes in LiCu202- We have used 
the projection procedure^ which is more accurate than 
the fitting procedure used in Ref. [H because the Wan- 
nier states in the former method are constructed from all 
electron DFT orbitals. The resulting Wannier orbitals 
of LiCu202, each centered at one Cu site, are shown in 
Fig. 5. One can see that the Wannier orbitals strongly 
overlap at oxygen and Cu + atoms. 

As we mentioned above, one of the possible microscopic 
mechanisms of exchange interaction is Hubbard-like AF 
superexchange. This interaction comes from hopping 
processes. We have calculated the hopping integrals be- 
tween orbitals of x 2 -y 2 symmetry in the Wannier function 
basis (Table I). The corresponding interaction paths are 



TABLE I: Calculated values of hopping parameters t%j be- 
tween x 2 — y 2 orbitals of copper atoms for the one-orbital 
model and estimated exchange interaction parameters using 
Eq.(22) (in meV). 
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presented in Fig. 6. One can see that for the largest hop- 
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FIG. 4: (Color online) Partial density of states obtained from 
LDA calculations. The blue solid and the dashed lines are the 
density of copper 3d states of x 2 — y 2 symmetry and the total 
density of 3d states, respectively. The green solid line is the 
density of oxygen 2p states. 



ping parameters we have good agreement with previous 
band fitting results. 12 However, there are also interaction 
paths which were not considered before. 

For nearest neighbors along the y axis, there is another 
contribution to the total exchange interaction, namely a 
FM "direct Coulomb exchange" between Wannier func- 
tions. The simplest physical representation of this sit- 
uation is presented in Fig. 7. In the case of the AFM 
configuration, the magnetic moment of the oxygen atom 
located between two copper atoms vanishes. By contrast, 
the magnetic moment of the oxygen atom in the FM case 
is not zero, and the energy gain is Jp x f3 4 , where j3 is the 



contribution of the oxygen atomic orbital to the Wannier 
orbital (see previous section). In order to calculate the 
couplings between magnetic moments, we use the follow- 
ing expression, which comes directly from the previous 
section: 



Ji 



<x*(U d - V dd ) + 



P 4 JpN 



(22) 



where N ox is the number of oxygen atoms on which 
the Wannier orbitals overlap. For simplicity, we ne- 
glect the intersite Coulomb interactions. The on-site 
Coulomb and intraatomic exchange interaction param- 
eters of copper atom are determined from first-principle 
constrained LDA calculations: U d =10 eV and =1 eV. 
Therefore the effective Coulomb interaction in Eq.(22) is 
U d = U d — = 9 eV. The value of the intraatomic ex- 
change interaction of oxygen atom, Jp , was estimated 
in LSDA+U calculations through the shift of oxygen 
2p band centers for spin- up, C\ and spin-down, C|: 
J* = (C T - C;)/M(0), where M{0) is the oxygen 
magnetization. The obtained value of 1.6 eV is in good 
agreement with previous estimations.— The value of a 2 
is related to the magnetization of copper atoms. Our 
LSDA+U results (see Table II) show that a 2 =0.58. Since 
the magnetic moment of the oxygen atom is the result of 
the magnetization of two copper atoms (see Fig. 7), /3 2 
= M(0)/2 = 0.09. Using Eq.(22) with the parameters 
defined above, one can calculate the exchange couplings 
between magnetic moments in LiCu202- These results 
are presented in Table I. One can see that the coupling 
between nearest neighbors along the y axis is strongly fer- 
romagnetic. These model considerations provide a micro- 
scopic explanation of the first-principle LSDA+U results 
presented in the next section. 



B. LiCu 2 2 : LSDA+U CALCULATION 




• O 

> Cu u 

f% en 2 




The analysis of the previous section shows that one 
should take into account Coulomb on-site correlations 




FIG. 5: (Color online) Wannier orbitals centered on neigh- 
boring copper atoms along y axis. 



FIG. 6: (Color online) The schematic representation of inter- 
action paths between copper atoms in LiCu2 02. 
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FIG. 7: (Color online) Schematic representation of AFM and 
FM configurations of magnetic moments at Wannier orbitals. 



and spin polarization of the oxygen atoms. This has 
been achieved using LSDA+U. The electronic structure 
of LiCii202 within LSDA+U is similar to that reported 
m Ref. 11. LiCu 2 2 is an insulator with an energy gap of 
0.7 eV. The values of the calculated magnetic moments 
(all values in units of are 0.58 for Cu 2+ and 0.18 for 
O. 

The next step of the investigation is a first-principle 
calculation of the isotropic exchange integrals for the 
Heisenberg model. In order to calculate the couplings 
between nearest neighbors along the y axis, one should 
use a method which takes into account the change of 
magnetization of the oxygen atoms. The most appro- 
priate one is the method-!- 9 in which the magnetic in- 
teraction Jij is estimated through the total energy dif- 
ference between the ferromagnetic and antiferromagnetic 
first-principle solutions (obtained, for instance, using the 
LSDA+U approach). The Heisenberg Hamiltonian de- 
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FIG. 8: (Color online) Ferromagnetic and antiferromagnetic 
configurations in an infinite chain with identical nearest neigh- 
bor exchange interaction J. The dashed line rectangle denotes 
the unit cell used for the mapping between the Heisenberg 
model and the first-principle LSDA+U approach. 

scribing the interaction between spins in the unit cell 
(Fig. 8) is given by 

H = 2zJSi ■ S 2 , (23) 
where z is number of nearest neighbors. The correspond- 



TABLE II: Results of LSDA+U calculations for the supercell 
1x2x1. M(Cu 2+ ) and M(O) are magnetic moments of 
copper and oxygen atoms, respectively. Etot is the relative 
total energy (in meV) normalized to 4 (number of copper 
pairs in unit cell). 



M(Cu^ + ) M(O) E total 



FM 
AFM 



0.58 
0.58 



0.18 
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ing total energies of the ferromagnetic and antiferromag- 
netic configurations of two spins are given by: 



E 



FM 



2zJS 2 



and 



E A 



FM 



= -2zJS 2 



(24) 



(25) 



Therefore the exchange interaction J is expressed in the 
following form: 



J 



E 



FM 



E 



AFM 



AzS 2 



(26) 



We have performed calculations for the (1x2x1) su- 
percell in ferromagnetic and antiferromagnetic configura- 
tions. The results are presented in Table II. One can see 
that in the case of the (1x2x1) supercell, the magneti- 
zation on oxygen atoms for the antiferromagnetic config- 
uration is zero, whereas the compensation of magnetiza- 
tion on oxygen atoms in the ferromagnetic configuration 
does not take place. Therefore first-principle LSDA+U 
calculations support the model considerations presented 
in the previous section. The ferromagnetic configuration 
has a lower energy than the antiferromagnetic one. Using 
Eq.(26) with z=2 and S=4, we get: JV. = -19.1 meV. 



TABLE III: Values of exchange interactions between mag- 
netic moments of LiCu2 02 compound (in meV). 
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In order to calculate the other magnetic interactions, 
we have implemented the Green's function method^ Ac- 
cording to this method, we determine the exchange inter- 
action parameter between copper atoms via the second 
variation of the total energy with respect to small devi- 
ations of the magnetic moments from the collinear mag- 
netic configuration. The exchange interaction parame- 



ters Jij for the Heisenberg model (Eq.(3)) with S=^ can 
be written in the following form 



20,21 



7T 



delm £ ( A i 



G 



m,m 
m ,m ' 
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where m (m', m" , m'") is the magnetic quantum number, 
the on-site potential A™" 1 ' = H™f - H™f and the 
Green's function is calculated in the following way 

„mn ju\ „m'n * ( h,\ 

G-'^=E \ En (27) 

k,n ' a 

Here cj™ n is a component of the nth eigenstate, and E™ 
is the corresponding eigenvalue. 

Our results are summarized in Table III. One can see 
that LSDA+U results are in good agreement with those 
obtained in our model analysis (see previous section) 
and disagree with previous theoretical estimates.— This 
agreement between the model analysis and LSDA+U re- 
sults is very encouraging, but clearly the ultimate test is 
to compare them with experiments. 

In that respect as well, the present results are a clear 
improvement with respect to previous estimates. Indeed, 
in contrast to results of paper Ref. [TJ], the ratio be- 
tween the strongest couplings J^/Jfj = -0.5 is in good 
agreement with the results of the neutron scattering ex- 
periments of Ref. IbH . This ratio is very important since 
it controls the pitch vector q of the helimagnetic state of 
LiCu202- The agreement is not perfect however: Our es- 
timates are about twice larger than the integrals deduced 
from experiments. Interestingly enough, it is possible, 
using the simple microscopic model Eq.(22), to identify 
the source of discrepancy between LSDA+U and exper- 
imental results. Indeed, the second term of Eq.(22) is 
very sensitive to the choice of 1 . For example, /3 2 =0.08 
leads to Jfj = -18.5 meV, which is in excellent agreement 
with LSDA+U results. For a 2 =0.65 and /3 2 =0.06, we 
obtain the following set of model exchange interactions: 
Jf, =-10 meV, Jff= 5 meV, J* =2.3 meV, J^=0.Q meV, 

J^ v =0A meV and J^ yz =0.5 meV. These model magnetic 
couplings are in good agreement with the experimental 
results. 

This proves the sensitivity of the results to the precise 
form of the Wannier functions. Now, it is well known that 
several sets of localized functions can be used to described 
a given band, 22 and the question of which Wannier func- 
tions should be used in the case of the determination of 
magnetic exchange, a point already raised by Anderson 
a long time ago (2 has not been settled yet. The present 



results call for further investigation of that issue. 

Another possible way however to improve the agree- 
ment between theory and experiment could be the follow- 
ing. From the experimental point of view, LiCu202 has 
spiral magnetic order in the ground state. Our study was 
performed for collinear magnetic configurations. There- 
fore, it would be more natural to calculate the exchange 
couplings using the magnetic structure observed exper- 
imentally. This goes beyond the scope of the present 
paper however. 

IV. DISCUSSION 

In conclusion, we have presented an analysis of the 
Hubbard model in the case of nearly 90° metal-oxygen- 
metal bonds dealing explicitly with Wannier orbitals. 
This has allowed us to derive an explicit expression of 
exchange integrals entirely in terms of parameters that 
can obtained from constrained LDA and LSDA+U. This 
expression can serve to interpret LSDA+U ab initio es- 
timates of the exchange integrals, and to establish their 
reliability. The analysis has been applied to the investi- 
gation of the magnetic couplings of LiCu202, allowing us 
to reach qualitative agreement with experiments, and to 
gain insight into the nature of exchange in that system. 
Because of the formation of a strongly hybridized, and 
energetically isolated combination of 3d x 2 _ y i and 2p or- 
bitals, a large moment is transferred to the O ions, and 
the magnetization of oxygen atoms has been proven to 
be the main source of ferromagnetism in LiCu202. 

We would like to thank I.V. Solovyev, A.I. Lichtcn- 
stein, A.O. Shorikov, J. Dorier and A. Gelle for helpful 
discussions. The hospitality of the Institute of Theoret- 
ical Physics of EPFL is gratefully acknowledged. This 
work is supported by INTAS Young Scientist Fellowship 
Program Ref. Nr. 04-83-3230, Netherlands Organiza- 
tion for Scientific Research through NWO 047.016.005, 
Russian Foundation for Basic Research grant RFFI 04- 
02-16096, RFFI 06-02-81017 and the grant program of 
President of Russian Federation Nr. MK-1573.2005.2. 
The calculations were performed on the computer cluster 
of "University Center of Parallel Computing" of USTU- 
UPI. We also acknowledge the financial support of the 
Swiss National Fund and of MaNEP. 



1 W.Heitler and F. London, Z. Physik 44, 455 (1927). 

2 W. Heisenberg, Z. Physik 49, 619 (1928). 

3 P.W. Anderson, Phys. Rev. 115, 2 (1959); Solid State 
Physics 14, 99 (Academic, New York 1963). 

4 J. Hubbard, Proc. R. Soc. London, Ser. A 276, 238 (1963) 

5 G.H. Wannier, Phys. Rev. 52, 191 (1937). 

6 J.J. Girerd, Y. Journaux and O. Kami, Chem. Phys. Lett. 
82, 534 (1981). 

7 C. Graaf, I. Moreira, F. Illas, O. Iglesias and A. Labarta, 
Phys. Rev. B 66, 014448 (2002). 



8 W. Ku, H. Rosner, W. E. Pickett, and R. T. Scalettar, 
Phys. Rev. Lett. 89, 167204 (2002). 

9 A. Huebsch, J. C. Lin, J. Pan, D. L. Cox, Phys. Rev. Lett. 
96, 196401 (2006). 

10 O. Gunnarsson, O.K. Andersen, O. Jepsen, and J. Zaanen, 
Phys. Rev. B 39, 1708 (1989). 

11 DA. Zatsepin el ai, Phys. Rev. B 57, 4377 (1998). 

12 A. A. Gippius et ai, Phys. Rev. B 70, 20406 (2004). 

13 T. Masuda et al, Phys. Rev. Lett. 92, 177201 (2004). 

14 T. Masuda et al, Phys. Rev. B 72, 014405 (2005). 



8 



O.K. Andersen, Z. Pawlowska and O. Jepsen, Phys. Rev. 
B 34, 5253 (1986). 

R. Berger, J. Alloys Compd. 184, 315 (1992). 

V.I. Anisimov et ai, Phys. Rev. B 71, 125119 (2005). 

I.I Mazin and D.J. Singh, Phys. Rev. B 56, 2556 (1997). 

Iberio de P.R. Moreira and F. Mas, Phys. Rev. B 55, 4129 

(1997). 

A.I. Lichtenstein, M.I. Katsnelson, V.P. Antropov, and 



V.A. Gubanov, J. Magn. Magn. Mater. 67, 65 (1987). 
V.V. Mazurenko and V.I. Anisimov, Phys. Rev. B 71, 
184434 (2005). 

See e.g. F. Lechermann, A. Georges, A. Poteryaev, S. Bier- 
mann, M. Posternak, A. Yamasaki, and O. K. Andersen, 
Phys. Rev. B 74, 125120 (2006), and references therein. 



